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We consider a two-dimensional dilute Bose gas above its superfluid transition temperature. We 
show that the t-matrix approximation corresponds to the leading set of diagrams in the dilute limit, 
provided the temperature is sufficiently larger than the superfluid transition temperature. Within 
this approximation, we give an explicit expression for the wave vector and frequency dependence 
of the self-energy, and calculate the corrections to the chemical potential and the effective mass 
arising from the interaction. We also argue that the breakdown of the t-matrix approximation, 
which occurs upon lowering the temperature, provides a simple criterion to estimate the superfluid 
critical temperature for the two-dimensional dilute Bose gas. The critical temperature identified by 
this criterion coincides with earlier results obtained by Popov and by Fisher and Hohenberg using 
different methods. Extension of this procedure to the three-dimensional case gives good agreement 
with recent Monte Carlo data. 
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PACS numbers: 05.30.Jp, 05.70.Fh, 03.75.Fi 

I. INTRODUCTION 

Renewed interest in two-dimensional (2D) superfluid 
systems has been recently prompted by the discovery of 
high-temperature superconductors. Before that, study- 
ing 2D bosonic systems was mainly motivated by ex- 
periments on adsorbed helium monolayers!] and spin- 
polarized hydrogen recombining on a helium film.B More 
recently, Bose condensation has been achieved experi- 
mentally in dilute gases of alkali atomscl and in atomic 
hydrogen,Q stimulating active research in this field both 
experimentally and theoretically.p 

Even though Bose-Einstein condensation is known not 
to occur at finite temperature for the ideal and inter- 
acting boson systems both in one and two dimensions, 
the absence of a condensate does not necessarily imply 
the lack of a phase transition tx> a superfluid state for 
an interacting 2D Bose system.B In this system, parti- 
cles with small momenta behave like a condensate and 
are responsible for the presence of a nonvanishing su- 
perfluid density.Q The critical temperature for the su- 
perfluid phase transition was estimated by Popovu us- 
ing a functional integral formalism. The same estimate 
for the critical temperature was later obtained by Fisher 
and Hohcnbergp using a renormalization-group approach, 
with the result 



27m 



mln In [l/(nrg)] 



(1) 



where n is the particle density, m the boson mass, and r 
the range of the interaction potential between bosons (wi 
set H — ks = 1 throughout). In these studies, Popovl 
and Fisher and HohenbergB approached the phase tran- 
sition from the superfluid phase and the normal phase, 
respectively. 



The present paper studies the 2D dilute Bose gas above 
its critical temperature T c by relying on conventional dia- 
grammatic methods. [The criterion for a 2D Bose system 
to be "dilute" will be specified below.] In this way, re- 
sults for thermodynamic quantities (like the critical tem- 
perature and the chemical potential) are most readily 
obtained. In addition, these results are amenable to ex- 
tension to more complex systems, such as the composite 
bosons occurring in the BCS to Bose-Einstein crossover 
problemB Finally, our method enables us to treat on 
equal footing the dilute Bose gas both in two and three 
dimensions. 

Previous studies of the 2D dilute Bose gas have ccpsick 
ered the superfluid phase-eiiher at zero temperaturet^rEJ 
or at finite temperature J3'£3 At finite temperature, the 
absence of a condensate in 2D required either the sep- 
aration of wave-vector integration into rapid and slow 
parts, or—the use of appropriate renormalization group 
methods.E-3 Previous approaches to the 2D dilute Bose 
gas, however, did not address the following issues: (i) 
The description of the normal phase (above the critical 
temperature) by standard diagrammatic methods, (ii) 
Upon lowering the temperature, the detection by these 
methods of a superfluid phase transition not associated 
with the establishing of long-range order. 

In this respect, standard many-body diagrammatic 
methods prove sufficient for a complete description of 
the normal phase of the 2D dilute Bose gas. In partic- 
ular, the t-matrix approximation for the self-energy will 
be shown to provide the correct description for a dilute 
Bose gas above T c , akin to the three-dimensional (3D) 
case. Both in 2D and in 3D, an explicit analytic expres- 
sion for the self-energy as a function of wave vector and 
frequency will be presented. 

In our approach, the expression (fil) for the critical 
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temperature obtained by PopovQ and by Fisher and 
Hohenbergu appears as a lower bound for the validity 
of the t-matrix as a controlled approximation for the di- 
lute Bose gas. The ensuing diagrammatic classification 
scheme for the dilute Bose gas will, in fact, be shown 
to break down when the temperature T approaches a 
lower temperature Tl, which coincides with the esti- 
mate (jl]) for the critical temperature given in Refs. 
and H In this way, the occurrence of a superfluid phase 
transition enters the diagrammatic theory for the di- 
lute Bose gas, since the physical mechanism leading to 
a breakdown of diagrammatic perturbation theory can 
only be the presence of a phase transition. In addi- 
tion, by applying our method to the 3D case we obtain 
{T L - T BB ) Abe ~ n l ^a, where T BE is the 3D Bose- 
Einstcin temperature and a is the scattering length., in 
agreement with recent Monte Carlo simulationstSEj on 
the 3D dilute Bose gas which yield for the critical tem- 
perature T c the same result we obtain for Tl- 

The paper is organised as follows. Section 2 sets up 
the diagrammatic theory for the 2D dilute Bose gas in 
the normal phase and calculates the corrections to the 
chemical potential and the effective mass due to the in- 
teraction. Section 3 discusses the occurrence of the su- 
perfluid phase transition for the dilute Bose gas through 
a breakdown of the approximations introduced in Sec- 
tion 2. The value of the breakdown-temperature Tl is 
then given both in two and three dimensions. Section 
4 gives our conclusions. In the Appendix, the full de- 
pendence of the t-matrix self-energy on wave vector and 
frequency is obtained both in 2D and in 3D by exploit- 
ing the diluteness condition. In addition, the validity of 
some approximations on which the theoretical arguments 
of the text rely is explicitly tested numerically. 

II. T-MATRIX APPROXIMATION FOR A 
DILUTE SYSTEM OF INTERACTING BOSONS 

In this Section, we analyze the diagrammatic theory 
for the 2D Bose gas in the normal phase and determine 
the leading contributions to the self-energy in the dilute 
limit. We give an analytic expression for the wave- vector 
and frequency dependence of the self-energy, which be- 
comes asymptotically exact in the dilute limit (as ver- 
ified in the Appendix). We further use the zero wave- 
vector and frequency value of this self-energy to dress 
the single-particle Green's functions of the theory in a 
self-consistent way. This step will enable us to identify 
a lower temperature Tl below which the classification of 
diagrams for the dilute Bose system breaks down, as dis- 
cussed in the next Section. In addition, we obtain the 
explicit leading corrections to the chemical potential fi 
due to_ interaction, thus recovering an earlier result by 
Popo\Jj, and to the effective mass. 



We begin by considering a 2D bosonic system inter- 
acting via a short-range two-body potential v(r) with a 
finite range ro, which becomes a S- function when ro — > 
(cf. Ref . |l2|) . We examine the dilute limit of this system, 
which is initially identified by the condition nr^ <C 1, 
where n is the bosonic density. [A stronger condition on 
the parameter nr\ will be required below.] We further 
consider temperatures above a nominal critical temper- 
ature T c but lower than an upper temperature of the 
order of 2™-= 2irn/m, at which quantum effects become 
important. til 

Under these assumptions, the selection of the diagrams 
yielding the leading contributions to the self-energy for 
the 2D dilute Bose gas in the normal state proceeds along 
similar lines as for the 3D dilute Bose gas. 

Akin to the 3D case, also in 2D every cycle (defined 
as a closed path constructed by a sequence of "bare" 
bosonic propagators, with a common wave vector and 
Matsubara frequency flow) introduces a Bose function 
( e (q /{ 2m )-n)/ T — 1) _1 , which appears after the summa- 
tion over the common frequency running along the cycle 
is performed. This function, in turn, cuts off the integral 
over the remaining wave-vector variable approximately 
at |q| ~ (mT) 1 / 2 , which is much smaller than the cut- 
off 1/t*o introduced by the potential (owing to the the 
diluteness condition and the assumption T ~ T n ). 

The description of the dilute Bose gas accordingly re- 
tains only those diagrams with a minimal number of cy- 
cles, that is, just one cycle. These diagrams are shown 
in Fig. 1(a) and constitute .tha so-called t-matrix approx- 
imation for the self-energy.l3iI3 




(b) (c) 

FIG. 1. (a) t-matrix approximation for the self-energy of 
an interacting Bose gas; (b) Particle-hole bubble; (c) Parti- 
cle-particle bubble. 

At T ~ T n higher-order diagrams can be estimated by 
replacing the bare potential by the t-matrix itself and as- 
signing to each additional cycle a factor ion p h(0). Here, 
to and n p h(0) are the zero frequency and wave- vector 
values of the t-matrix and of the particle-hole bubble 
(Fig. 1(b)), respectively. It will further shown below that 
the product ion p h(0) cx 1/ ln[l/(wo)] <C 1 in the dilute 
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limit. 

We will verify, however, that the classification of dia- 
grams based on the cycle argument breaks down upon 
lowering the temperature down to a value Tl- More 
precisely, we will find that the t-matrix correctly de- 
scribes the dilute Bose gas in the temperature range 
Tl <T n , in the sense that no other diagrams besides 
the t-matrix itself need to be included. 

The self-energy corresponding to the set of diagrams 
depicted in Fig. 1(a) reads: 
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where the t-matrix t(p, p',P) is defined by the integral 
equation 

t(p,p',P)=v(p-p')-Tj2 J p«-q) 

xG (r f ) G (i +? ) t(q,p,,p) (3) 

with the notation q = (q, iu> u ) (u) v = "^jll/T - v integer - 
being a bosonic Matsubara frequency) S3 

To lowest order in the density, all single-particle 
Green's functions G in Eqs. (Q) and (^) are considered 
to be bare ones. Quite generally, self-energy insertions 
in the Green's functions become relevant when the phase 
transition is approached. It is shown in the Appendix 
that, in the dilute limit, the dependence of the self-energy 
Sj(g) on q can be disregarded when calculating physical 
quantities, so that one may set S t (g) ~ S t (0) in all single- 
particle Green's functions entering Eqs. (||) and (|3|) and 
re- absorb the constant E t (0) by a shift of the chemical 
potential. In this way, when the chemical potential is 
expressed in terms of n and T, one obtains 



fx - S t (0) = no 



(4) 



where fio = HoIm, T) is the chemical potential of the 2D 
ideal Bose gas.EJ 

In the following, we shall consider all Green's functions 
in Eqs. (||) and (|J) to be self-consistently dressed by the 
self-energy (|^) with q = 0, resulting in an "improved" t- 
matrix approximation which can be adopted to approach 
the critical temperature more closely. 

The explicit expression of E t (q) for arbitrary values of 
q is given by Eq. (A7) of the Appendix, which is asymp- 
totically valid in the dilute limit. From that expression, 
the shift of the chemical potential due to interaction as 
well as the relevant effective mass can be obtained. For 
the chemical potential, it is sufficient tq_know the value 
of E t (0), which from Eq. flA7|) becomes :El 
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By entering in this exoression the analytic form of the 
chemical potential /io we obtain eventually 



St(0) 



iim/m 



ln(l/[m7>2|ln(l-e- T "/ T )|]) 
87rn/m 



ln[l/(nr§)] + ln[T n /(2vrT)] - In | ln(l - e~^/ T )\ 
^ 8tto 

~ mln[l/(nrg)] l J 

the last result holding for TJ ln[l/(nrg)] < T < T ni 
which includes the temperature range (Tl ^ T < T n ) 
of physical interest. Equation @ provides the leading 
self-energy term for the 2D dilute Bose gas in the normal 
state, from which the shift of the chemical potential is 
obtained as ji = fiQ + S t (0). 

Note that the expression @ is temperature indepen- 
dent in the temperature range Tl ^ T < T n we are 
considering. The only temperature dependence of /i thus 
originates from fiQ. In particular, at T = T c we substi- 
tute the value of T c from Eq. (0) onto the expression Mm, 
yielding: 



fi(T c ) = ^ (T C ) + E t (0) 
2im 
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provided that lnln(l/(rtrQ) is sufficiently larger than 
unity (dilute limit). This result coincides with the value 
of the chemical potential obtained in Refs. ^,|| where the 
critical temperature was approached from below. 

The effective mass can eventually be calculated after 
analytic continuation of St (q). The result is 
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rn 



41nln[l/(nrg)] 



(8) 



Details of the derivation of Eq. 
Appendix. 



are reported in the 



III. BREAKDOWN OF THE BOSONIC 
T-MATRIX APPROXIMATION AND THE 
SUPERFLUID PHASE TRANSITION 

We pass now to show that the selection of diagrams 
made in the previous Section, which was based on the 
cycle argument, breaks down upon lowering the tempera- 
ture when the superfluid phase transition is approached. 
Specifically, consideration of the temperature at which 
the particle-hole diagrams (which where discarded by the 
cycle argument) are no longer negligible in comparison 



3 



with the particle-particle diagrams, will lead us to iden- 
tify a lower temperature Tl that turns out to coincide 
with the critical temperature ([[]) determined in Refs. |?j 
and ||. 

To determine the range of validity of the t-matrix ap- 
proximation, it is enough to compare the particle-particle 
bubble Upp of Fig. 1(c) (which constitutes the building 
block of the t-matrix of Fig. 1(a)) with the particle- 
hole bubble Ilph of Fig. 1(b). This statement follows 
from the classification of diagrams we have made because 
t ~ npp(O)^ 1 in the dilute limit [cf. Eqs. and (JA4|)], 
yielding n p h(0)/n pp (0) as the small parameter of the the- 
ory at T »j T n . The ratio n p h(0)/n pp (0) grows, however, 
upon lowering the temperature below T n and reaches the 
value of unity at the temperature Tl introduced above. 

The particle-particle bubble is given by 



n PP (o) = T£ p 



d 2 q 

(2tt) 2 



G(q) G(-q) 



(9) 



where G(q) — (iuj u — q 2 /(2m) + /i ) 1 according to the 



arguments above. With the notation £(q) 
jLto, we obtain: 



|q| 2 /(2m)- 
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where ng(i) = (e x > — 1) is the Bose function and the 
last asymptotic equality holds in the dilute limit nr 2 <C 1 
and for temperatures of the order of T n (such that \/j,q\ ~ 

T n ). 

The particle-hole bubble for q — is, as usual, given by 
n p h(0) = dn/dfi. For the two dimensional Bose gas the 
chemical potential Ik known analytically for all tempera- 
tures and densitiesEJ, such that 



dn 



2tt 1 - e»°/ T 



(11) 



At T ~ T n , /i (n, T n ) ~ T n and dn/d[i ~ m, such 
that the ratio n ph (0)/n p p(0) w 1/ ln[l/(nr 2 )] < 1 in 
the dilute limit. When T <C T n , on the other hand, 
/io| <C T and dn/d/j, ~ mT/(2ir\fj,o\), which is much 
larger than the corresponding value at T ~ T n . A lower 
temperature Tl can thus be reached, such that the ratio 
n p h(0)/n pp (0) equals unity when 
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Entering now the asymptotic expression |jL*o | ~ Te T ™/ T , 
which holds for T -C T„, yields eventually 
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T n 



lnln[l/(w 2 )] 



(13) 



valid under the assumption lnln[l/(nr 2 )] ^> 1 (which 
defines the diluteness condition in 2D). This expres- 
sion for Tl coincides with the estimate ([[]) for the crit- 
ical temperature given by PopovQ and by Fisher and 
Hohenbergu. Note that the double-log dependence of Tl 
on nr§ originates, on the one hand, from the log depen- 
dence of n pp (0) on \/j,q\ and, on the other hand, from 
the exponential dependence of fiQ on T at low tempera- 
tures. Note also that the more stringent diluteness con- 
dition lnln[l/(nr 2 ) ;§> 1 (in the place of the original 
nr 2 <C 1) is required to get a finite temperature range 
(Tl J$ T < Tn), where self-energy diagrams can be se- 
lected by the diluteness condition. 

On physical grounds, the only mechanism which may 
lead to the breakdown of the diagrammatic classification 
in terms of a small parameter (as explicitly seen above) 
is the occurrence of a phase transition. The breakdown 
temperature Tl is thus expected to provide an estimate 
for the superfluid critical temperature T c , also because 
the two temperatures Tj, and T c are expected to have 
the same functional dependence on nr 2 , as comparison 
of our expression for Tl with the value (jl]) for the critical 
temperature obtained in Refs. |?],|| indeed shows. 

Our identification of Tl with T c is further confirmed 
by applying the same sort of arguments to the 3D diluie. 
Bose gas, for which Monte Carlo results are available Jl3't£l 
In this case, the particle-particle and particle-hole bub- 
bles are given, respectively, by 



n PP (o) 



m 



2ir 2 ro 

where a is the scattering length, and 

2mn 1 / 3 T BE 



n P „<»> = | 



3C(3/2) 4 /3 t _ T BE 



(14) 



(15) 



since n ^ -9C(3/2) 2 (T - TBE) 2 /(167rT B E) for the chem- 
ical potential of the ideal Bose gas (valid near the Bose- 
Einstein temperature T BE = 27r/C(3/2) 2 / 3 n 2 / 3 /m). The 
temperature at which these contributions coincide then 
defines the 3D breakdown temperature Tf° . One ob- 
tains: 
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3C(3/2) 4 / 3 



2.33 n 



V3, 



(16) 



This result agrees with recent Monte Carlo simulations 
for the 3D hard-core Bose gas in the dilute limit ,t£l which 
yielded precisely (T C -T BE )/T BE = (2.3±0.25)n 1 / 3 a, and 
coincides with the analytic result of Ref. [22], which was 
obtained by a completely different method. 

The result ( |T6| ) should be regarded as altogether non 
trivial, since previous analytic treatments of the 3D di- 
lute Bose gas resulted either in different dependences of 
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(T c — Tbe)/7be on the parameter n 1 / 3 a, e.g., of the type 
(n^a) 1 / 2 [cf. Ref. || or (n^a) 2 / 3 [cf. Ref. M, or in 
the same linear dependence on the parameter ?iL| a, but 
with a different proportionality coefficient .E3cZl In our 
approach, the linear dependence on the parameter n 1 / 3 a 
of the temperature shift has been directly related to the 
quadratic dependence of the free-boson chemical poten- 
tial on T — Tbe near Tbe in 3D. 



IV. CONCLUDING REMARKS 

In this paper, we have considered the two-dimensional 
dilute Bose gas in the normal phase, in the interesting 
temperature region ranging from an upper temperature 
T n (below which quantum effects become important) to a 
lower temperature Tl (which we have identified as the su- 
perfluid critical temperature). In this temperature region 
we have analyzed the ordinary diagrammatic theory and 
organized it in powers of the parameter 1/ fn[l/(wo)], 
which was assumed to be small compared to unity. 

In this way, the standard t-matrix has been identified 
as yielding the dominant set of diagrams for the self- 
energy when nr 2 <C 1. Further analysis of the theory 
to define the temperature range where the t-matrix ap- 
proximation holds has, however, led us to consider the 
stronger condition lnln[l/(na 2 )] ^ 1 as characteristic of 
the "dilute" Bose gas in two dimensions, thus confimiing 
the criterion introduced by Fisher and Hohenberga via 
different methods. 

Our identification of the lower temperature Tjr, (and, 
thus, of the superfluid critical temperature) rests on the 
finding that the diagrammatic classification scheme for 
the 2D dilute Bose gas breaks down at this lower temper- 
ature, in the sense that additional diagrams (besides the 
t-matrix) become also important at Tl and the hierarchy 
established for the dilute gas no longer holds. In this re- 
spect, it may be worth mentioning that our criterion to 
identify the critical temperature does not contradict the 
usual criterion which defines T c (n) as the temperature 
where the equation [x(n, T c ) = £(g = 0; n, T c ) is satisfied. 
Solving, in fact, for this equation requires one to rely 
on an approximation for the self-energy which is valid 
even at T c . We have seen, however, that in our case the 
t-matrix approximation for the self-energy breaks down 
before reaching T c , as soon as the critical region above T c 
is approached. 

The finding that an estimate of the superfluid critical 
temperature can be obtained from the ordinary diagram- 
matic theory in the normal phase, both in two and three 
dimensions, constitutes per se a nontrivial result, espe- 
cially because the nature of the (superfluid) transition 
in two and three dimensions is quite different (involving, 
respectively, quasi-long-range order and true long-range 
order). In addition, our approach is rather straightfor- 
ward and amenable to direct implementation to more 



complex physical situations. 

In this respect, a possible application of our results 
may be the normal state of high-temperature cuprate 
superconductors, which are quasi-two-dimensional sys- 
tems. Experiments related to the normalc§E2l and 
superconducting^ state in these systems suggest, in 
fact, that a correct description of their properties might 
require an intermediate (crossover) approach between 
the Fermi liquid theory (weak-coupling) and the di- 
lute Bose gas approach (strong-coupling). For this rea- 
son, crossover theories have been considered by sev- 
eral authors, both for the normal stateEH and the 
broken-symmetry phase.HEj Since a reliable study of the 
crossover problem should require a good knowledge of at 
least the extreme (weak- and strong-coupling) limits, the 
approach developed in this paper mighLshed light on the 
strong-coupling limit of these theories.EI 
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APPENDIX A: 

In this Appendix, we obtain an analytic expression for 
the t-matrix self-energy, which is valid above Tl in the 
dilute limit. From this expression, the effective mass (at 
q = 0) is calculated. In addition, the approximation 



E t (g) ~ E t (0) 



(Al) 



(which is crucial for the arguments presented in the text) 
will be justified theoretically and checked numerically. 
Both two-dimensional and three-dimensional cases will 
be considered for completeness. 



1. Two dimensions 

It is convenient to parametrize the short-range poten- 
tial by a separable potential in wave-vector space, by 
setting u(k — k') = VQWkWk> with — 6(ko — |k|) and 
k = r^ 1 . In this way, Eq. (|J) can be readily solved to 
yield i(p,p',P) = w p w p 't(P), with 



t -. ( P)=»o- + TE/|^ 



(A2) 



where G(q) — — q 2 /2m + Ho) 1 , consistently with 
the assumption ( Al) and the ensuing Eq. (0). This choice 
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of G(q) will lead us to verify the key assumption ( |AlD in 
a self-consistent manner. The frequency sum in Eg. ( |A2| ) 
can be performed explicitly, yielding 



(2^) 2 



P 



-q)G 



d 2 q 2 1 + n B (C P /2- q ) +n-s(Cp/ 2 +q) 



(2tt) 2 q P 2 /(4m) + q 2 /m-2/i -ift,. 



(A3) 



In the dilute limit (nr 2 < 1) and for T L < T < T n , 
the Bose functions appearing in the above expression can 
be neglected, since they yield contributions smaller by a 
factor 2>g <C 1 with respect to the term retained. The 
integration over the wave vector q then yields 



1 

Vq 
m 

— In 

47T 



k 2 /m + 2|/xq| + P 2 /(4to) - ifl v 
2\fj. \ +P 2 /(4m) -iQ u 



(A4) 



where In stands for the principal branch of the complex 
logarithm. The t-matrix self-energy is obtained by insert- 
ing expression (A4) into Eq. (g) which, for a separable 
potential, becomes 



E t (q,^) = -2T 

x G(q',a, 



E 



d 2 q' 2 
(2tt) 2 W (i'-i)/ 2 

w„)f(q + q', SV 



(A5) 



To perform the frequency sum in Eq. ( A5) w e ex ploit 
the analytic properties of t(P,z). From Eq. (A4), af- 
ter the replacement i£l v > — > z, it can be readily verified 
that i(P, z) has a simple pole for z — P 2 /(4m) + 2|/io| + 
k 2 /[m(l — e" 1 /" )] (with ?)o = rnuo/(47r)) and a branch 
cut along the real axis for P 2 /(4m) + 2|^ | < Re(z) < 
P 2 /(4m) + 2|/xo| + A: 2 /TO. The frequency sum in Eq. ( |A5| ) 
can be then performed by a contour integration, yielding 
three distinct contributions: one from the simple pole of 
the Green's function G(q', z), one from the simple pole of 
the t-matrix, and one from the integration along the cut 
of f (q + q', z). The term originating from the simple pole 
of the t-matrix, which occurs for Re(z) > k 2 /m, is expo- 
nentially suppressed by the Bose factor l/(exp(/3z) — 1) 
and is thus negligible. The term associated with the inte- 
gration along the cut can be estimated to be smaller than 
the term from the simple pole of the Green's function by 
a factor l/ln[l/(nr 2 )], and is also negligible in the dilute 
limit. We are thus left with the expression: 
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E t (q,w„)=2 
fc 2 /m 



d 2 p 



V2 



Mo I 
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- q 2 /(2m 



p 2 /(4m) 



1 m 
— + T- 

Vq 47T 
n -1 



I A*o | + q 2 /(2?7i) - p 2 /(4m) - iuj v 



the range of p over which the log-term in the denomina- 
tor of Eq. (A6) varies appreciably. When performing the 
integration over p in Eq. ( |A6| ) , we can then approximate 
p = q in the logarithm, which is in this way factored out 
of the integral, yielding the following asymptotic expres- 
sion for the self-energy: 



S«(g) 
2n 



V2 



m 

— In 
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1 

VQ 



kl/m + |/i | + q 2 /(4m) - iu> v 



|Mo| + q 2 /(4m) - iu v 



(A7) 



The analytic expression (|A7| ) can be compared with 
the numerical calculation for the t-matrix self-energy, ob- 
tained by retaining all contributions to Eq. (A5). 




(a) 




(A6) 



(b) 

FIG. 2. OJu) (units of T n ) for a 2D dilute Bose gas at 

T — T„ and nrl = 1.6 x 10~ 3 . Numerical results for v = 10 
(diamonds) and v — 1 00 (squares) are compared with the 



Note that the Bose factor in Eq. (A6) is peaked about 
p = q with a width of order T 1 / 2 , which is smaller than 



analytic expression (A7) (full and dotted lines, respectively). 



We have found that the asymptotic expression (A7) 
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reproduces extremely well the numerical results for ( |A5| ) 
when the diluteness parameter is sufficiently small. As 
an example, in Fig. 2 the analytic expression for 
is compared with the numerical results for the choice of 
parameters T = T n , mv /(Air) — 1, and nr 2 — 1.6x 10~ 3 . 
One can see that in this case the agreement is excellent. 



The expression (A7) for the self-energy can also be 
used to verify that, when evaluating physical quanti- 
ties (such as the density n), the wave- vector and fre- 
quency dependence of S t (g) is actually irrelevant, so 
that one can approximate £ t (g) — S t (0) as anticipated 
in Eq. (Al). Thi s is because the presence of the loga- 



rithm in Eq. (A7) makes the dependence of S t (q) on q 
rather slow. The approximation T, t (q) ~ S 4 (0) is thus 
justified over a large portion of q space and can be ex- 
ploited to evaluate physical quantities. In particular, for 
T = T n we have obtained numerically that the relative 
error when evaluating n = —T J (2^1 G{q) alterna- 



tively with G(q) = [iQ v — q 2 /(2m) - 



■M-St(g)]- 1 



with 



G(q) -v G (q) = [ift v - q 2 /(2m) + n - Sf(O)]- 1 is less 
than 1% for nrl < 10" 2 . 



Finally, the t- matrix self-energy (A7) can be exploited 
to calculate the effective mass for the dilute Bose gas. 
Recall that, once the retarded selfpepergy is known, the 
effective mass can be calculated asE3 



i9ReE(|q|,aO 
dui 



9ReS(|q|,w) 



9|q| 



(A8) 



where the derivatives are meant to be calculated at the 
quasi-particle pole defined by the equation 

~2 



/i + RcS(q,w q ) 

2m 



(A9) 



In general, the effective mass m* depends on q. Here we 
are interested in its value at q = 0, which is relevant at 
low temperatures. 



The self-energy (A7) can be analytically continued via 
the replac eme nt iu> n — > u) + i0 + . The quasi-particle-pole 
equation (A9) at q = can then be solved asymptoti- 
cally, to yield 

,lnln[l/(nr 2 )r 



uj = - M +S t (0) 1-2 



ln(l/(nr 2 )] 



(A10) 



as it can be verified by inserting the value (AlC) for coo in 
the quasi-particle-pole equation and by discarding sub- 
leading terms in the dilute limit nr\ <C 1. The deriva- 
tives of S(q, u>) at (q = 0,ui = ojq) can also be readily 
calculated. The effective mass at q = is then given by 

m* _ l + l/(21nln[l/(rarg)]) 
m 



1 - l/(41nln[l/(w 2 )]) 



l + l/(41nln[l/(nr 2 )]). 



(All) 



Note the occurrence of the same double-log dependence 
characteristic of the temperature Tl ■ 



2. Three dimensions 

The three-dimensional case can be treated in a parallel 
fashion to the two-dimensional case. By considering the 
same separable potential adopted in 2D, and by following 
the same steps which lead to Eq. (A4), we now obtain: 



vq 2tt z 



x arctan ■ 



,1/2 



-2^o 



P 2 

iil v 

4m 



ko/m 



1/2 



'^0 + ^ 



(A12) 



where the complex arctan is defined, as usual, in terms of 
the principal branch of the complex logaritm as follows: 



1 1 — iz 
arctan z = — m 

2 1 + iz 



(A13) 



Like in 2D, the t-matrix t(P, z) has a branch cut along 
the real axis for P 2 /(4m) + 2\fi \ < Re(z) < P 2 /(4m) + 
2|// |+&o/ m > an d a simple pole located along the real axis 
for Kez > P 2 /(4m) + 2|^ | + k\jm. Upon trasforming 
the frequency sum in ( A5) into a contour integration, the 
t-matrix simple pole contribution will be again strongly 
suppressed by the Bose factor, while the term associated 
with the integration along the cut can now be proven to 
be smaller than the term originating from the simple pole 
of the Green's function by a factor (nr-p) 1 / 2 . In the dilute 
limit, only the simple pole of the Green's function there- 
fore contributes to the contour integration. By the same 
argument leading to Eq. (|A7j), we obtain the following 
asymptotic expression for the t-matrix self-energy: 



St(g) 
2n 



V2 



,.V2 



1 

arctan 



2tt 2 



fco 



Mo I 



4m 
-1 



1/2 



Mo 



4m 



(A14) 



The asymptotic expression ( A14 ) h as a lso been checked 
against numerical calculation of Eq. ( A5 ) in three dimen- 
sions. In Fig. 3 the analytic expression ( A14 ) is compared 
with the numerical results for the choice of parameters 
T = T n , (27r 2 )/(mw fco) -> 0, and n 1 / 3 ^ = 1 x 1CT 2 . 
Even in 3D the agreement is excellent. 

The approximation £((<?) — S(0) has further been 
checked by evaluating the particle density. In this case, 
we have found that the error introduced by the approx- 
imation (Al) in the estimate for the density is less than 
1% when nrl $i 5 x 10~ 3 . 

Finally, the solution of the quasi-particle-pole equa- 
tion (A9) at q = in 3D is given by 



^0 = -M + s t(0) = |mo| , 



(A15) 
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while the retarded self-energy near the quasi-particle pole 
is given by 



E(q,w) = 



2)i 



7T , i i 7r/(2aA(q,^)) + l 
2 2 111 7r/(2aA(q,w))-l 



(A16) 



with A(q, a;) = yjuj — q 2 / (4m) — |/io| and where we 
have introduced the scattering length a via the relation 
m/(47ra) f» (mk ) / (2tt 2 ) . By using the definition (AS), 
the effective mass at q = can be readily calculated, 
leading to the result: 



to* Sirna 3 ( 4 

to to \ ir z 

Anna 3 ( 4 
wl + 1 + — 

TO V 7T^ 



47rna 3 / 4 
1 + 1 + — 

TO \ 7T Z 



(A17) 



Note that the interaction increases the (quasi)-particle 
mass both in 3D and in 2D with respect to its bare value. 




(a) 



-0,05 - 



:-0,l - 



-0,15 - 




-0,2 



e-a p 



[]■■■□■ a- • -b 



Q..Q---QQ- 



0,5 



1,5 

q/k () 

(b) 



2.5 



FIG. 3. Y,(q,u u ) (units of T n ) for a 3D dilute Bose gas at 
T = T n and n 1/3 r = 1 x 10~ 2 . Numerical results for v = 10 
(diamonds) and v — 10 (squares) are compared with the 
analytic expression ( A.14) (full and dotted lines, respectively). 
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